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I. INTRODUCTION 


The physics of exotic nuclei with unusual N/Z ratios (isospin) has attracted world-wide 
attention. These nuclei are usually loosely bound and often exhibit resonances with a pro- 
nounced single-particle character since the Fermi surface approaches the continuum. In some 
nuclei, the valence nucleons can be easily scattered into single-particle resonant states (or 
Gamow states) in the continuum which is essential for the existence of exotic phenomena. 
For example, the halo in Li [1] can be successfully reproduced [2] and giant halo nuclei [3] 
are predicted by the relativistic continuum Hartree-Bogoliubov (RCHB) theory [4] in which 
the contribution from the continuum has been proven to be crucial to understand these halo 
phenomena. There are also other attempts to include the contribution from the continuum 
in the non-relativistic [5] and relativistic approaches [6]. By taking into account a few reso- 
nant states close to the continuum threshold with the BCS method [6], the relativistic mean 
field (RMF) theory seems to be able to reproduce the RCHB results and describe the main 
effects of the continuum. It has also been shown that the contribution of the continuum to 
the giant resonances mainly comes from Gamow states [7, 8]. Therefore, a proper treatment 
of resonant states is important for a deeper understanding of the properties in exotic nuclei. 

Presently, several techniques have been developed to study resonant states in the con- 
tinuum. One of them is the R-matrix theory in which resonance parameters (i.e. energy 
and width) can be reasonably determined from fitting the available experimental data [9]. 
The extended R-matrix theory [10] and the K-matrix theory [11] have also been developed. 
The conventional scattering theory is also an efficient tool for studying resonances. More 
precisely, the phase shift method is commonly used to determine parameters for a resonant 
state which corresponds to a pole of the S matrix [12]. By discretizing the continuum, 
the contribution of the resonant states can be self-consistently taken into account via a 
Bogoliubov transformation in coordinate space [2, 13). 

Computationally, it is desired to deduce the properties of unbound states from the eigen- 
values and eigenfunctions of Hamiltonians for bound states so that the methods developed 
for bound states can still be used. For this purpose, the bound-state-type methods have 
been developed, including the real stabilization method [14], the complex scaling method 
(CSM) [15], and the analytic continuation in the coupling constant (ACCC) method [16]. 


In the real stabilization method, the continuum is discretized in a box and the resonance 


can be found according to the fact that its energy should be stable against the change of 
the box size. Sometimes it is difficult to find the appropriate box sizes, and it also turns out 
that this approach is suitable for narrow resonances only [17]. In the CSM, the wave func- 
tion for a resonant state is transformed to the square-integrable wave function for a bound 
state by rotating the position vector in coordinate space into the complex plane. Since the 
CSM involves the calculations of complex matrix elements, the corresponding difficulties 
related to complex eigenvalue problem still exist. In the ACCC method, a resonant state 
becomes a bound state if one increases the attractive potential and the energy, width and 
wave function for the resonant state can be obtained by an analytic continuation carried 
out via a Padé approximant (PA) from the bound-state solutions. Compared with other 
bound-state-type methods, the ACCC approach is very effective and numerically quite sim- 
ple because many methods available for bound-state problems can be used and the PA for 
analytic continuation can be implemented easily. 

Within the framework of the non-relativistic Schrodinger equation, the ACCC method 
has been applied to investigate the energies and widths for resonant states in light nuclei 
combined with few-body methods [17, 18], and to study single-particle resonant states in 
spherical and deformed nuclei by solving the Schrodinger equation with Woods-Saxon po- 
tentials [19]. Based on the Schrödinger and Dirac equations, the stability and convergence 
of the energies and widths for single-particle resonant states with square-well, harmonic- 
oscillator and Woods-Saxon potentials have been investigated and their dependence on the 
coupling constant interval and the order of the PA have been discussed in detail [20, 21]. 

The relativistic mean field (RMF) theory is widely used to describe properties of nuclear 
matter and finite nuclei successfully, not only for those nuclei near the valley of stability, 
but also for exotic nuclei with large neutron (proton) excess [2, 22]. Previously the ACCC 
method was employed to study energies and widths for resonant states in stable nuclei 160 
and Ca in RMF theory [23]. In the present work, we aim at exploring the single-particle 
resonant states, including resonance parameters and wave functions, in exotic nuclei by the 
ACCC method within the framework of the RMF theory. This paper represents the first 
application of the RMF-ACCC method to obtain both the wave functions and resonance 
parameters. The paper is organized as follows. In Sec. II the ACCC method combined with 
the RMF theory is presented. The numerical details are given in Sec. HI and its application 


to the doubly magic nucleus '?Zr in Sec. IV. Finally, we give a brief summary in Sec. V. 


Il. THEORETICAL FRAMEWORK 
A. Relativistic Mean Field Theory 


The basic ansatz of the RMF theory is a Lagrangian density whereby nucleons are de- 
scribed as Dirac particles which interact via the exchange of various mesons (the scalar o, 


vector w and iso-vector vector p) and the photon [4] 
L=wv(id-—M)p + $0,000 —U(o)—- OO” 
+3m2w,w" — Ru R” + imp, p” — FF F” (1) 
—~gswow =T guppy T IpbPTw > ep Ay, 


where M is the nucleon mass and m, (go), Mw (gu), and Mm, (gp) are the masses 


(coupling constants) of the respective mesons. A nonlinear scalar self-interaction 
1 
U(o) = goo” + Lo + Fo of the ø meson has been included [24]. The field 


tensors for the vector mesons are given as 
QY = olw” — owt, 
R” = Hp’ — a" pl — gP(p" x p”), (2) 
Fey = OPA” — OV A". 


The classical variation principle gives the following equations of motion 


lo-p+VW(r) + B(M + Vs(r)) bi = exis (3) 


for the nucleon spinors, where c; and Y; are the single-particle energy and spinor wave 


function respectively, and the Klein-Gordon equations 
(—Ao T U'(c)) = JoPs, 
(A + mZ)w? = gaj"(r), 


(=A + m?) p” — gt’ (r), 


—A Ap(r) ejst), 


for the mesons, where 


Volt) = gush + gT + Zell = 73), 


V(r) 


(5) 


IJo%,; 


are the vector and scalar potentials respectively and the source terms for the mesons are 


A _ 
Ps S 2 Vitir 
3 i 
Hr) = Debi hi, 
Fo (6) 
- O = Do WYT 
A. {uo 
HO) = Vd 


It should be noted that the contribution of negative energy states are neglected, i.e., the 


vacuum is not polarized. Moreover, the mean field approximation is carried out via replacing 


meson field operators in Eq. (3) by their expectation values, since the coupled equations 
Eq. (3) and Eq. (4) are nonlinear quantum field equations and their exact solutions are very 
complicated. In this way, the nucleons are assumed to move independently in the classical 
= meson fields. The coupled equations can be solved self-consistently by iteration. 

ee For spherical nuclei, the potential of the nucleon and the sources of meson fields depend 
= only on the radial coordinate r. The spinor is characterized by the angular momentum 
quantum numbers l, j, m, the isospin t = +i for neutron and proton respectively, and the 


other quantum number i. The Dirac spinor has the form 


xXa(t), (7) 


where Yla (0, @) are the spinor spherical harmonics. The radial equation of the spinor, i.e. 


Eq. (3), can be reduced to [4] 


(CS + ÉRC) + (M+ VG E) = G(r), 
(È + EGY) — [M — Vnt) FË) = eM), (8) 


in which V,(r) = Ww (r) + Vs(r), Vin(r) = Ww (r) — Vs(r), and « = (—1))t+4+/2(5 4.1/2). The 


meson field equations can be reduced to 


o? 20 
CE td) O k 


where ¢ = g, w, p, and photon ( mg = 0 for photon). The source terms read 


and 


Arr? pelr) 


—goPs — 920° (r) — gs0? (r) for the o field, 


JwPo for the w field, 
(10) 

IpP3(r) for the p field, 
ep-(r) for the Coulomb field, 

` 2 2 
= 2 (G:(r)| — |Filr)|?), 

A 2 2 
= DGr) + Flr), 

z (11) 
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n=] 


IMs 


Ms 


(IGp(r)/? + |Fp(r)I?)- 


3 
ll 
mn 


B. Analytic Continuation in the Coupling Constant for Resonance Parameters 


and Gamow Wave Functions 


By solving Eqs. (8) and (9) in a meshed box of size Ro self-consistently, one can calculate 


the ground state properties of a nucleus. The vector potential V(r) and the scalar potential 


Vs(r), energies and wave functions for bound states are also obtained. By increasing the 


attractive potential as V,(r) — AV,(r), a resonant state will be lowered and becomes a 


bound state if the coupling constant A is large enough. Near the branch point Ao, defined 
by the scattering threshold k(\o) = 0 [16], the wave number k(A) behaves as 


iy ÀA — ào, 1>0, 
K(X) ~ (12) 


W= Jai 


These properties suggest an analytic continuation of the wave number k in the complex A 

plane from the bound-state region into the resonance region by Padé approximant of the 

second kind (PAII) [16] 

cot ee + eon? +...+ CLT? 

kla) x pe) = r l l 
(7) (7) 1+ diz + dəz? +... + dna’ 


(13) 
where z = VA— ào, and co, c1,...,CL, d1, d2,...,dy are the coefficients of PA. These co- 
efficients can be determined by a set of reference points v; and k(x;) obtained from the 
Dirac equation with A; > Ao, i = 1,2,..,0+ N + 1. With the complex wave number 
k(\ = 1) = k, + tk;, the resonance energy Æ and the width F can be extracted from the 
T 
relation £ = E — is (E,T € R) and k? = e? — M?, i.e., 
(M2 + k2 — k?)? + 4k2k? + (M° + k? — k?) 
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T = 4/2y/(M2 + k2 — k2)2 + 4k2k2 — 2(M2 + k2 — k2). (14) 


In the non-relativistic limit (k « M), Eq. (14) reduces to 


ke k? 


E r- 2krki 


2M ’ M ` 


(15) 


It is evident that the continuation in the coupling constant can be replaced by the con- 
tinuation in k plane along the k(A) trajectory determined by Eq. (13) to the point kpr 
corresponding to the wave number for Gamow state, i.e., kg = k#NI(A = 1) [16]. Similarly, 
the wave function y(kp,r) for a resonant state can be obtained by an analytic continuation 
of the bound-state wave function y(k;,r) in the complex k plane. One can also prove that 
the wave function y(k,r) is an analytic function of the wave number k in the inner region 
r < Ro where the Jost function analyticity dominates [16]. Therefore, we use the technique, 
which has been adopted to find the complex resonance energy, to determine the resonance 
wave function y(kr,r). Firstly, we construct the PA to define the resonance wave function 


at any point r in the inner region (r < Ro) [16] 


NI(p p) — Pil(k,r) 2 ag(r) +ai(r)k + a(r)k? +... +ar(r)k? 
P (k, j= Qulk,r)  1+b(r)k+ba(r)k?+...+by(r)kð ’ = 


where the coefficients a;(r) (i = 0,1,..., L) and b;(r) (j = 1,2,..., N) are dependent on r. 


These coefficients can be determined by a set of reference points k; and ọ(k;,r) obtained 


from the Dirac equation with A; > Ap, (i = 1,2,...,L+N-+1). The resonance wave function 
lkr, r) = pl (kp, r) (r < Ro) can be extrapolated in this way. 

Eq. (8) can be rewritten as two decoupled Schrédinger-like equations for the upper and 
the lower components respectively [25]. For neutrons, the Schrédinger-like equation for the 
upper component reads 

POCO) leid ahte) + (e? MAGE) <0, a7 
in the outer region where V(r) ~ 0 and Vs(r) ~ 0. Its solution is the well known Riccati- 


Hankel function 


A*(z) = Alz) E ijkl), z= kr, (18) 


K 


where j,(z) = zj,(z) is the usual regular solution, i.e. Riccati-Bessel function, and fi,(z) = 
zn,(z) is the irregular solution, i.e. Riccati-Neumann function. The outer wave function is 


matched to the inner wave function at r = rm < Ro 


A 


p(k tm) = O (kn) [Ju(krerm) + D(ke)ina(krerm)| , (19) 


where C(kp) is the coefficient for matching and D(kr) = tand,(kr) with (kr) the phase 
shift. It has been found that 6,(kg) is almost a constant when rm is large enough. Given 
the upper component, the lower component F),(r) can be calculated from the relationship 


between the upper and the lower components derived from Eq. (8), i.e., 


Pale) = FSP [ira lenr) + Din)in-alkar)]- (20) 


Finally, the resonance wave function is normalized according to the Zel’dovich proce- 


dures [16]. 


II. NUMERICAL DETAILS 


In this section we give the details on how the ACCC method is realized within the frame- 
work of the RMF theory and exploit the analyticity of the eigenvalues and eigenfunctions of 
the Dirac equation with respect to the coupling constant. In RCHB calculations °°Ca is the 
doorway for exotic phenomena such as giant halos [3], in which the neutron resonant states 
play an important role in the continuum. We select the neutron resonant state 71g9/2 in °°Ca 


to investigate the analyticity of the eigenvalue and eigenfunction for the Dirac equation with 
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respect to the coupling constant. The RMF equations are solved in a spherical box of the size 
Ro = 20 fm, with a mesh size of ôr = 0.05 fm [26]. Several well used effective interactions 
for the RMF Lagrangian will be used and the corresponding results will be compared. The 
energy, width and wave function for the v1gg9/2 resonant state in °°Ca are obtained with the 
procedure given in the last section. As the stability and convergence (e.g., the dependence 
on PA order) of the energies and widths in the Dirac and Schrédinger equations have been 
investigated in Refs. [20, 21], we do not discuss them here and use L = N = 4 PA for the 
present calculations. 

In Tab. I, we present the energies and widths for v1g9/2 in Ca calculated with effective 
interactions NL3 [27], NLSH [28], TMA [29] and the newly developed PK1 [30] parametriza- 
tion. The results from the scattering phase shift method (RMF-S) obtained similarly as in 
Ref. [6] are also included for comparison. The first step in RMF-S is to solve the Dirac 
equation [Eq. (8)] with scattering boundary conditions for the continuum spectrum. The 
solutions, upper and lower components, have the same form as here in Eq. (19) and Eq. (20) 
respectively. The coefficients C and D in Eq. (19) and Eq. (20) are fixed by the normaliza- 
tion condition of the scattering wave functions and the phase shift ô is calculated from the 
matching conditions. It should be mentioned that the energy used in the wave functions is 
real and not unique, as opposed to complex and the unique value which is determined by the 
ACCC method in advance. In the vicinity of an isolated resonance, the resonance energy 


and width are extracted from the derivative of the phase shift ô in Breit-Wigner form, i.e., 


d5(E) T/2 
dE (E, -B+E E 


After discretizing the real energy for scattering solutions and phase shift variation, the 
resonance energy is uniquely determined. 

One can find from Tab. I that the results from the ACCC and the scattering methods 
agree well. Different effective interactions give similar energies and widths. Therefore in the 
following calculations, we use only the effective interaction NL3 for illustration. 

The resonance energy E and width I for the neutron resonant state v1gg/2 in ®°Ca, as 
a function of the coupling constant A, are exhibited in Fig. 1. Filled circles and crosses are 
the solutions of the Dirac equation [Eq. (8)], and the former is used as input in the ACCC 
method. Solid curves are the outcome of L = N = 4 PA [Eq. (13)]. The dashed vertical 
line corresponds to the branch point Ay) > 1. When A > Xo, the particle is bound. With 


A decreasing from Ag down to 1, a nonvanishing imaginary part appears in the complex 
energy expressed by Eq. (14). At \ = 1, one gets the resonance energy E and width Tr. The 
agreement between the exact solutions of the Dirac equation with large coupling constant 
(crosses in Fig. 1) and the extrapolated results of PA (solid energy curve) is quite good, 
which demonstrates the validity of the ACCC method. The smooth solid line indicates good 
analyticity of the eigenvalue for the Dirac equation with respect to the coupling constant. 
The same results can also be obtained for the case of L = N = 5 and L = N = 6 PAs, 
similar as in Ref. [20, 21]. 

In Fig. 2, we show the real and imaginary parts of the upper component G(r) and the 
lower component F'(r) at r = 5 fm and r = 10 fm respectively as functions of the resonance 
energy E (which corresponds to the coupling constant A displayed in Fig. 1). Symbols in 
Fig. 2 have the same meaning as those in Fig. 1. The neutron feels much stronger effect from 
the nuclear potential at r = 5 fm than at r = 10 fm. Similar to the case of the eigenvalue 
discussed above, the real and the imaginary part of each component behave smoothly, which 
indicates good analyticity of the eigenfunction for the Dirac equation with respect to the 
coupling constant. 

The wave function obtained from the analytic continuation is connected continuously to 
the free wave function at the matching point r,,. On one hand, the matching point r,, 
should be large enough to make sure that the potential at rm vanishes, which is required 
by the outer part of the wave function. On the other hand, rm should be far enough from 
the box radius Ro to guarantee that the inner part of the wave function can be obtained 
accurately. In the lower panel of Fig. 3, the potential V(r) + Vs(r) for °°Ca is plotted. One 
finds that when r > 14 fm, the potential becomes smaller than 1074 MeV which is good 
enough for the precision in our calculations. The upper limit of the matching point rm can 
be estimated by the behavior of the density from the upper panel in Fig. 3. The neutron 
density pn(r) in Ca and Prig, (r) = [|G(r)|? + |F(r)|?]/4rr? decline rapidly in the vicinity 
of 20 fm because of the boundary condition in the RMF calculation. The density curves 
in the range r < 18 fm are smooth and behave reasonably. Therefore, the matching point 
should be chosen at any point between 14 fm and 18 fm. We present the wave function 
for the neutron resonant state v1g9/2 in °°Ca with different matching points in Fig. 4. As 
expected, the wave function is the same when changing the matching point rm from 14 fm 


to 18 fm. In the following calculations, the wave function will be matched at rm = 18 fm, 
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unless otherwise specified. 


IV. APPLICATION TO SINGLE-PARTICLE RESONANT STATES IN !”°ZR 


We apply the ACCC method within the framework of the RMF theory to calculate the 
energies, widths and wave functions for the resonant states in 1??Zr. Since !?*Zr is the core of 
the giant neutron halo predicted in RCHB calculations [3], and the resonant states in '°?Zr 
have been obtained by the scattering method [6], we choose this nucleus for comparison. 

In Fig. 5, the energies and widths for the resonant states v3p3/2, v3pi/2, V2f7/2, V2fs/2, 
vlhgj2 and 1213/2 are presented as a planar E-P plot. The results of the ACCC and the 
scattering methods are in good agreement with each other for most of these states. From 
both methods, v2fs/2 and 1113/2 have large widths, while v2 f7/2 and v1hg/2 are very narrow. 
The width for v3p3/2 from the scattering method is slightly larger than our calculation. As 
for the resonant state v3pj;/2, neither the energy nor the width can be extracted from the 
scattering calculation. From a quantum mechanical point of view, these single-particle 
resonant states are quasi-stationary ones captured by centrifugal barriers. The decay width 
for a resonant state can be roughly explained by the penetration through the barrier. For 
those states with the same J, i.e., the same centrifugal barrier, the higher state has a larger 
width, e.g., for the two f states v2f5/2 and v2 f7/2. Although v1hg/2 is higher than v2 fs/2, 
its width is smaller because of higher centrifugal barrier. A similar argument also holds for 
a much higher but narrow state 1113/2. 

The wave functions for the resonant states v3p3/2, V2 f7/2, V2fs/2 and v1hg/2 are respec- 
tively exhibited in Figs. 6, 7, 8, and 9. Since the imaginary part of the wave function is 
neglected in the scattering calculations [5, 6], only the real parts of the upper and the lower 
components are given and compared in these figures. In Fig. 6, one finds a considerable 
difference for the upper components at large r in v3p3/2 between these two methods. This 
difference is consistent with the difference in energy and width as seen in Fig. 5. A similar 
difference is also seen for v2f5/2 in Fig. 8. For those states where the two methods give 
nearly the same energies and widths, e.g., v2 f7/2 and v1hg/2, the wave functions agree quite 
well. The small difference for the wave functions between these two methods found in Figs. 7 
and 9 might be from the neglect of the imaginary part of the wave function in the scattering 


method. 
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V. SUMMARY 


In summary, the method of the analytic continuation in the coupling constant (ACCC) 
has been developed within the framework of the relativistic mean field (RMF) theory and 
used to investigate the resonance energies, widths and wave functions for single-particle 
resonant states in exotic nuclei. The analyticity of these quantities, as functions of the 
coupling constant, are exploited for the neutron resonant state v1gg/2 in Ca. The energies, 
widths and wave functions for some resonant states in '2*Zr are calculated and compared 
with those obtained from the scattering method. The results from the ACCC method agrees 
satisfactorily with those from the scattering calculation as well as available data [23]. It is 
expected that the RMF-ACCC method, combined with the BCS theory, can be generalized 


to describe exotic phenomena in unstable nuclei. 
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Tab. I: Energies and widths (E, T) in MeV for the neutron resonant state v1gg/2 in Ca calculated 


by the ACCC (RMF-ACCC) and the scattering (RMF-S) method with the effective interactions 
NL3, NLSH, PK1, and TMA. 


NL3 NLSH PK1 TMA 


RMF-ACCC [MeV]| (0.49, 0.0005) (0.47, 0.0004) (0.77, 0.0009) (0.77, 0.0009) 


RME-S [MeV] 


(0.54, 0.0001) (0.53, 0.0001) (0.85, 0.0006) (0.83, 0.0005) 


Energy and width [MeV] 


E 1.5 1.4 H 1.1 1.0 0.9 


Fig. 1: Energy E and width I for the neutron resonant state v1gg/2 in 60Ca, as a function of 
the coupling constant A. Filled circles and crosses are the solutions of the Dirac equation, and the 


former is used as input in the ACCC method. Solid curves are the outcome of L = N = 4 PA. 
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Fig. 2: The real and imaginary part of the upper component G(r) and the lower component F'(r) 
for the neutron v1gg/2 state in °°Ca at (a) r = 5 fm and (b) r = 10 fm as functions of the resonance 
energy E (which corresponds to the coupling constant À displayed in Fig. 1). Symbols have the 


same meaning as those in Fig. 1. 
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Fig. 3: Neutron densities for Ca and for the neutron v1gg j2 State (upper panel); the potential 
W(r) + Vs(r), the single-particle neutron bound states and the neutron resonant state v1gg/2 


(lower panel). 
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Fig. 4: Real parts of the upper and lower components of the radial wave function for the neutron 


v1gg/2 state in 60Ca calculated by the ACCC method with different matching points rm. 
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Fig. 5: Energies and widths for the neutron state v3p3/2, v3pi/2, V2f7/2, V2f5/2, Vihg/2, and 
v1tj3/2 in 1222r. Solid circles represent the results of the ACCC method, while open circles denote 


the results of the scattering method. 
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Fig. 6: Real parts of the upper and lower components of radial wave functions for the neutron 
resonant state v3p3/2 in 12277. Solid and dashed curves represent the results of the ACCC and the 


scattering method respectively. 
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Fig. 7: Similar as Fig. 6, but for the neutron resonant state v2f7/2 in 1227r. 
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Fig. 8: Similar as Fig. 6, but for the neutron resonant state v2f5/2 in ee 
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Fig. 9: Similar as Fig. 6, but for the neutron resonant state v1hg/2 in 1227r, 
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